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ABSTRACT 

>• ' Wide-field, very long baseline interferometry (VLBI) observations of the 

^ ■ nearby starburst galaxy NGC 4945, obtained with the Australian Long Baseline 



in 

O ■ Array (LEA), have produced 2.3 GHz images over two epochs with a maximum 

angular resolution of 15 mas (0.3 pc). 15 sources were detected, 13 of which corre- 
spond to sources identified in higher frequency (3 cm and 12 mm) ATCA images. 
Four of the sources are resolved into shell-like structures ranging between 60 and 
110 mas (1.1 to 2.1 pc) in diameter. From these data the spectra of 13 compact 

■ radio sources in NGC 4945 were modelled; nine were found to be consistent with 

■ free-free absorbed power laws and four with a simple power law spectrum. The 
free-free opacity is highest toward the nucleus but varies significantly through- 
out the nuclear region (tq ~ 6 — 23), implying that the overall structure of the 
ionised medium is clumpy. Of the 13 sources, 10 have steep intrinsic spectra 
associated with synchrotron emission from supernova remnants, the remaining 
sources have fiat intrinsic spectra which may be associated with thermal radio 
emission. A non-thermal source with a jet-like morphology is detected ~ 1" 
from the assumed location of the AGN. A type II supernova rate upper limit 
of 15.3 yr~^ is determined for the inner 250 pc region of the galaxy at the 95% 
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confidence level, based on the lack of detection of new sources in observations 
spanning 1.9 years and a simple model for the evolution of supernova remnants. 
A type II supernova rate of > 0.1(f/10^) yr~^ is implied from estimates of su- 
pernova remnant source counts, sizes and expansion rates, where v is the radial 
expansion velocity of the supernova remnant in km s~^. A star formation rate 
of 2.4(t;/10^) < SFR{M > SMq) < 370 yr'^ has been estimated directly 
from the supernova rate limits and is of the same order of magnitude as rates 
determined from integrated FIR (1.5 Mq yr~^) and radio luminosities (14.4 ± 1.4 
((5/8.8) M0 yr^^). The supernova rates and star formation rates determined for 
NGC 4945 are comparable to those of NGC 253 and M82. 

Subject headings: galaxies: active - galaxies: individual (NGC 4945) - galaxies: 
starburst - radiation mechanisms:general - supernovae: general - techniques: 
interferometric 



Introduction 



This paper is the second of a series of papers on the sub-parsec scale properties of local 
{D < 10 Mpc), bright (51.4 > 10 mJy) southern (5 < —20°) starburst galaxies. Starburst 
galaxies are the result of processes acting on a wide ra nge of scales and the mselves drive 
energetic activity on a wide range of scales. In Paper I (ILenc fc Tingayl l2006l ) of this series 
we discussed the properties of the nuclear starburst region of NGC 253. Direct observations 
of supernova remnants enabled an investigation of the star-formation and supernova history. 
Furthermore, high-resolution radio observations were used to map the ionized gaseous envi- 
ronment of the starburst region. These observations provide a link between the large-scale 
dynamical effects that feed the starburst region, the activity in the starburst region, and the 
energetic phenomenon that in turn is driven by the starburst. In the present work (Paper 
II of the series) we discuss the sub-parsec scale properties of NGC 4945 and implications on 
the star-formation and supernova rate of this galaxy. 

NGC 4945 is a nearly edge-on sp iral galaxy associate d with the Centaurus group of 
galaxies at a declination of ~ —40° ( IWebster et al.Ml979l ). The ga laxy has been classi- 



fied as SB(s)cd or SAB(s)cd and has an optical extent of 17 arcmin (Ide Vaucouleurs 



fied as SB(s)cd or 


Braatz et al. 


1997) 


1996 


) and 8.1 Mpc 



1964 



the magnitude of stars at the tip of the red giant branch using Hubble Space Telescope (HST) 
Advanced Camera for Surveys (ACS) images and arrived at a more reliable distance estimate 
of 3.82 ± 0.31 Mpc. This new distance estimate will be adopted throughout this paper and 
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implies a spatial scale of ~ 19 pc arcsec 



While it is classified as a Seyfert 2 galaxy (IBraatz et al.l 119971: iMadeiski et al.l l2000l ) 



it exhibits both starburst and Seyfert characteristics (e.g. IWhiteoakI Il986l ). The galaxy 
is one of the strongest and ric hest sources of extragalatic molecular fines with ~ 50 line 



features detected to date (e.g. iHenkel et al.l Il994j : ICurran et al.l l200ll : IWang et al.l 12004 ). 



Evidence of an active galacti c nucleus (AGN) comes from strong and variable fiard x-ray 
emission (llwasawa et al.lll993r) . H2O megamaser emission fias been observed from tfie nucleus 
( iGreenfiill et al.l 119971 ) and is assumed to be distributed in a ~ 3 pc disk around tfie AGN. 



Altfioug h tfie central region o f NGC 4945 is fiigfily obscured by gas and dust at optical 
wavelengtfis (jMadejski et al. 2000 ) , it is one of the three brightest IRAS point sources beyond 
the Magellanic Clouds (jSanders et al" 2003). Neutral hydrogen observations of NGC 4945 
reveal the si gnature of a ba r associated with the disc and can be traced out to a radius 
of ~ 7 kpc Jott et al.ll200ll l C0(2— 1) o bservations sugges t that the bar extends inwards 



to within ~ 100 pc of the galaxy centre (|C 



and dust into fast-rotating disk of cold gas (iDahlem et al.l 



lou et al.ll2007f). where it may act to feed ga s 



I993I: iMauersberger et al. 



HNC(l-O) maps of the disk reveal a dense inner region (jCunningham fc WhiteoaklboOS) that 



may fuel a circumnuclear starburst ri n g of s imilar extent observed by iMarconi et al.l (2000) 



996|). 



with HST PAa images. IChou et al.l (120071 ) also detect a smaller-scale (r < 27 pc), high 



density, kinetically de-coupled component in the nu clear region and suggest that it may be 
evidence of a circumnuclear torus around the AGN (lAntonuccil Il993l ) . 



Associated with the nuclear region, and having its origin in the inner starburst region, 
is a powerful wind that forms a cone-shaped plume perpendicular to the disk of the galaxy. 
The p lume, extending ~ 500 p c north-west, is clearly observed i n Ha spectral line o bserva- 



tions (IRossa fc Dettmaii 120031 ) and Chandra soft x-ray images ( Scfiurcfi et al.l 12002). Sucfi 



Strickland et al 



plumes fiave been found in otfier nearby starburst galaxies (e.g. 

are believed to be driven by supernovae in tfie starburst region (jPoane &: Matfiew^ 
Strickland! [200! . 



2004) and 



1993 



Tfie bar associated witfi tfie disk of tfie galaxy, a fast-rotating disk of cold gas, a starburst 
ring and tfie conical nuclear outflow are all cfiaracteristic s tfiat are observed in otfier nearb y 
starburst galaxies sucfi as NGC 253 (ILenc fc Tingayll2006l ) and M82 (IForster Scfireibeiil2000l ). 
Tfiese features are interpreted in terms of bar- induced gas dynamics and star formation. Even 
tfiougfi NGC 4945, NGC 253 and M82 are associated witfi different galaxy groups and are 
situated in different locations in tfie celestial spfiere, tfiey sfiare several furtfier similarities. 
All tfiree galaxies are nearly edge-on, all fiave distances close to 4 Mpc and all fiave similar 
infrared luminosities (iRice et al.lll988l ). It is tfiese similarities tfiat make tfiem particularly 
interesting targets for study and comparison in otfier respects. 
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Figured] presents images of NGC 4945 at various wavelengths, on different spatial scales, 
to illustrate the relationships between the distribution of cool stars (K-band), ionised hydro- 
gen gas (Ha), the inner starburst region (radio and Paa), and the wind eminating from that 
region (Ha and x-ray). The supernova remnants that are a primary subject of this paper 
are concentrated within the inner 250 pc of the galaxy. 

The observation of new supernovae and their associated remnants in the starburst re- 
gions of nearby galaxies is rare as these regions are highly obscured by gas and dust. Only 
a few have been observed in local starbursts and are generally revealed in x-rays where such 
events are more conspicuous and then followed up with multi-wavelength observations. In a 
few cases these can be traced back to an earher optical or infrared outburst, the presumed 
superr iova event, as with SN 1996cr in Circ inus (IBauer et al.l 120071 ) and SN 1978K in NGC 
1313 (jSmith et al.l 120071 : iRyder et al.lll993l ). However, in most cases, where the remnants 
are embedded deep in the starburst region and only x-ray and radio observations can reveal 
their presence, they are discovered long after the event and their ages remain uncertain . 



Such a population of supernova remn ants has been discovered in M82 (jPedlar et al.lll999 



McDonald et al.ll2002l ) and NGC 253 flTingavl I2004J : iLenc fc Tingavl l2006l ) using VLBI ob- 
servations. 

Motivated by the existence of compact radio sources and strong indicators of free-free 
absorption in similar starburst galaxies, such as NGC 253 and M82, we have embarked on 
a program of VLBI observations of NGC 4945 and other prominant Southern Hemisphere 
starburst galaxies. This paper reports results from new observations of the nuclear region 
of NGC 4945 at a range of frequencies between 2.3 GHz and 23 GHz. The resulting data 
have better resolution and sensitivity than previous observations and we have used them to 
investigate in detail the free-free absorbed spectra of supernova remnants and H II regions 

O) 



in the NGC 4945 starburst (§ 14. ip . In particular we are interested in probing the structure 
of the ionised environment of the supernova remnants, constraining the supernova rate in 
NGC 4945 using our new ATCA and VLBI images (§ l4.5p and re-evaluating estimates of the 
star formation rate based on our results (§ 14. 6p . 

A subsidiary focus of this project is to achieve high sensitivity, high fidelity, high reso- 
lution, and computationally efficient imaging over wide fields of view (at least compared to 
traditional VLBI imaging techniques). As such, the associated exploration of the relevant 
techniques is a small first step toward the much larger task of imaging with the next gener - 
ation of large radio telesc opes, for example the Square Kilometre Array (SKA: iHalll (120051 ): 
Carilli fc Rawlingj (120041 )). The SKA will use baselines ranging up to approximately 3000 
km and cover fields of view greater than one square degree instantaneously at frequencies 
of around 1 GHz. The computational load required for this type of imaging task is vast 
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and will drive the development of novel imaging and calibration algorithms and the use of 
super-computing facilities. Investigating the performance of these techniques now, under the 
most challenging current observing conditions, is therefore a useful activity. 



2. Observations, correlation and data reduction 

2.1. LB A Observations 

A VLBI observation of NGC 4945 was made on 13/14 May, 2005 using a number of the 
Long Baseline Array (LEA) telescopes: the 70 m NASA Deep Space Network (DSN) antenna 
at Tidbinbilla; the 64 m antenna of the Australia Telescope National Facility (ATNF) near 
Parkes; 5 x 22 m antennas of the ATNF Australia Telescope Compact Array (ATCA) near 
Narrabri used as a phased array; the ATNF Mopra 22 m antenna near Coonabarabran; the 
University of Tasmania's 26 m antenna near Hobart; and the University of Tasmania's 30 



m an tenna near Ceduna. The observation utilised the S2 recording system (jCannon et al. 



19971 ) to record 2 x 16 MHz bands (digitally filtered 2-bit samples) in the frequency ranges: 
2252 - 2268 MHz and 2268 - 2284 MHz. Both bands were upper side band and right circular 
polarisation. 

A second epoch LBA observation was made on 21 March, 2007 using the same set of 
radio telescopes except that only 3 x 22 m antennas of the ATCA were used as a phased 
array to maintain a wide field of view given the longer baseline configuration of the ATCA 
at that time. The second epoch observation utilised hard disk data recorders (Phillips et al. 
2008, in preparation); allowing dual circular polarisation data across 4 x 16 MHz bands to be 
recorded at Parkes, ATCA, and Mopra, and 2 x 16 MHz bands at the remaining antennas. 
Observing parameters associated with each of the LBA observations are shown in Table [H 

During each of the VLBI observations, three minute scans of NGC 4945 {a = 13'^5™27!50; 
S = — 49°28'6'.'00 [J2000]) were scheduled, alternating with three minute scans of a nearby 
phase reference calibration source, J1237-5046 (a = 12^^37"^15!239286; 6 = -50°46'23'.'17260 
[J2000]). 



2.2. LBA correlation 



The first epoch data were correlate d using the ATNF Long Baseline Array (LBA) pro- 
cessor at ATNF headquarters in Sydney (IWilson et al.lll992l ). The data were correlated using 
an integration time of 2 seconds and with 32 frequency channels across each 16 MHz band 
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(channel widths of 0.5 MHz) to hmit the effect s of bandwidth smearing (a form of chro- 
matic aberration) and time- averaging smearing (jCottonlll999l : iBridle fc Schwablll999l ). The 
ATCA primary beam hmits the field of view of this observation to a half-width half-maximum 
(HWHM) of ~ 45". At the HWHM point, bandwidth smearing and time-averaging smearing 
losses are estimated to be approximately 15% and 4%, respectively. 



The second epoch data were correlated using the DiFX Software Correlator (jPeller et al. 



20071 ). The data were correlated using an integration time of 2 seconds and with 64 frequency 
channels across each 16 MHz band (channel widths of 0.25 MHz). The ATCA primary beam 
limits the field of view of this observation to a HWHM of ~ 1'. At the HWHM point, 
bandwidth smearing and time-averaging smearing losses are estimated to be approximately 
8% and 7%, respectively. 



2.3. LBA Data Reduction 

The correlated data were imported into the AIP^ package for initial processing. The 
data for the phase reference source were fringe-fit (AIPS task FRING) using a one minute 
solution interval, finding independent solutions for each of the 16 MHz bands. The delay 
and phase solutions for the phase reference source were examined and averaged over each 
three minute calibrator scan, following editing of bad solutions, before being applied to both 
the phase reference source and NGC 4945. Further flagging of the data was undertaken via 
application of a flag file that reflected the times during which each of the antennas were 
expected to be slewing, or time ranges that contained known bad data. Finally, data from 
the first 30 seconds of each scan from baselines involving the ATCA or Parkes were flagged, 
to eliminate known corruption of the data at the start of each scan at these two telescopes. 

During correlation, nominal (constant) system temperatures (in Jy) for each antenna 
were applied to the correlation coefficients in order to roughly calibrate the visibility am- 
plitudes (mainly to ensure roughly correct weights during fringe-fitting). Following fringe- 
fitting, the nominal calibration was refined by collecting and applying the antenna system 
temperatures (in K) measured during the observation, along with the most recently measured 
gain (in Jy/K) for each antenna. Further refinement in the calibration was complicated by the 
complex structure of the phase reference source J1237— 5046. To account for this structure a 



^The Astronomical Image Processing System (AIPS) was developed and is maintained by the National 
Radio Astronomy Observatory, which is operated by Associated Universities, Inc., under co-operative agree- 
ment with the National Science Foundation 
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DIFMAP ( Shepherd. Pearson and Taylor 1994 ) task, cordum^ f Lenc et aL 2006 ). was 



new 

developed to enable the transfer of all phase and amplitude corrections made in DIFMAP 
during the imaging process to an AIPS compatible solution table. The phase reference data 
were averaged in frequency and exported to DIFMAP where several iterations of modelling 
and self- calibration of both phases and amplitudes were performed. The resulting contour 
maps of the phase reference source are shown in Figure [2] and have a measured one sigma 
RMS noise of 0.342 and 0.292 mJy beam~^ in epochs 1 and 2, respectively, cordump was 
then used to transfer the phase and amplitude corrections back to the unaveraged AIPS data 
set. 

Averaging in frequency was performed in AIPS to increase the channel widths to 2 
MHz; this significantly reduced the data set to a manageable size and improved imaging 
performance. While this also reduced the overall field of view, the starburst region of interest 
was unaffected by the averaging process. To facilitate imaging in DIFMAP, a ParselTongu^ 
script was written to convert the averaged frequency channels into intermediate frequencies 
(IFs). This conversion allowed DIFMAP to treat the frequency channels independently in 
the (u,v) plane rather than averaging them together, thus avoiding any further bandwidth 
smearing effects during the imaging process. 

Initial imaging of NGC 4945 was performed at reduced resolution (80 x 32 mas) by ex- 
cluding data from the Hobart and Ceduna antennas (the longest baselines); uniform weight- 
ing was applied to maximise resolution at the expense of image noise. Several iterations of 
CLEAN and phase self-calibration were performed in DIFMAP. The resulting images are 
shown in Figure [3] and achieve a one sigma RMS noise of 0.152 and 0.082 mJy beam~^ for 
epoch 1 and 2, respectively. The image map parameters for both figures are shown in Table 

m 

A second phase of imaging was repeated at full resolution (16 x 15 mas) with data 
from all antennas. The imaging process followed the same procedure to that used in the 
initial imaging phase except that natural weighting was applied to provide a compromise 
between the increased resolution available with the full array and the available image noise. 
Image maps for both epochs were created for each of the sources that fell above the six sigma 
detection limit. The images from both epochs were then combined to produce a high dynamic 



^The cordump patch is available for DIFMAP at|http: //astronomy . swin. edu. au/$\sim$elenc/DifinapPatches/ 

•^A Python scripting tool for AIPS. ParselTongue was developed in the context of the ALBUS 
project, which has benefited from research funding from the European Community's sixth Frame- 
work Programme under RadioNet R113CT 2003 5058187. ParselTongue is available for download at 
[http : / / www . radionet- eu . org/ rnwiki/PaxselTongue 
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range image map for each of the detected sources. The epoch 1, epoch 2 and combined image 
maps for each of the detected sources are shown in Figure [6] and the associated map statistics 
tabulated in Table [21 The epoch 1, epoch 2 and combined maps, achieve a one sigma RMS 
noise of 0.149, 0.110 and 0.075 mJy beam~^, respectively. These figures compare favourably 
with the expected theoretically predicted thermal noise of ~ 0.120, ~ 0.080 and ~ 0.070 
mJy beam~^, respectivelj0. 



2.4. ATCA Observations and Data Reduction 

Observations of NGC 4945 were carried out using all six antennas in the EW367 and 6C 
configurations of the Australia Telescope Compact Array (ATCA) on the 2006 March 18, 24 
and 25. The EW367 configuration is compact with most baselines less than 367 m in length, 
the 6C configuration is for high resolution observations and has a spread of baselines up to 
6 km in length. The flux scale was set using PKS B1934— 638 and J1326— 5256 was used 
as a phase calibrator. During each observing run, we alternated between two setups every 
6.5 minutes. The first setup observed at 17 GHz and 19 GHz and the second setup at 21 
GHz and 23 GHz. For intermediate frequencies of 4.8 GHz and 8.64 GHz, archive data was 
obtained from the Australia Telescope Online Archive. The data were observed with the 6A 
configuration of the ATCA, which has a spread of baselines up to 6 km in length, on 1993 
November 19 and 20. The flux scale was set using PKS B1934-638 and PKS 1320-446 was 
used as a phase calibrator. A list of all ATCA observations reported in this paper and their 
associated observing parameters are tabulated in Table [TJ 



All ATCA data were initially calibrated using the MIRIAD software package (jSault. Teuben. fc Wright 



1995h. Subsequent calibration, decon volution and imaging was performed with the DIFMAP 



( jshepherd. Pearson and Taylor 1994) software package and final images were created using 



the KARMA software package (iGoochI Il996l ). The resulting contour maps are shown in 



Figure [5] with all associated beam parameters, image rms noise, peak and integrated flux 
densities tabulated in Table [2J 



''Estimated with the ATNF VLBI sensitivity calculator: ^http : //www . atnf . csiro . au/vlbi/calculator/ 
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3. Results 

3.1. Identification of Sources 

A list of sources detected above the 6a detection threshold in both epochs of the short- 
baseline VLBl data is given in Table [31 The 6a threshold provides a conservative false 
detection rate of less than one in 3000. Flux density errors of ±10% are list ed due to 



uncertainties in the absolute flux density scale for Southern Hemisphere VLBI (IReynolds 



19941 ) ■ The total flux density of each source was determined by summing its CLEAN model 
components in the image plane. The sources were also modelled with elliptical Gaussian 
components to provide an approximate measure of the full-width half-maximum (FWHM) 
size, the fltted size and position angle (PA) of each source is also listed in Table [31 Two 
of the detected sources, 27.5291 — 04.631 and 27.6164 — 03.373, exhibit structure that is 
signiflcantly larger than the beam width and are shown in greater detail in Figure [H Source 
27.5291 — 04.631, in particular, appears to exhibit a jet-like morphology and has been divided 
into three components in Table [31 27.5291 - 04.631, 27.5402 - 04.720 and 27.5607-04.779. 

Based on the reduced resolution images, a total of 12 discrete sources fell above the 6a 
threshold (0.91 mJy beam~^) of the flrst epoch observation of NGC 4945, whereas 13 fell 
above the 6a threshold (0.51 mJy beam~^) of the second epoch observation. When imaged 
with the longer baselines included (Hobart and Ceduna), eight of these sources were detected 
in the flrst epoch data and nine in the second epoch. While source 27.4646—05.174 did not 
achieve the 6a detection threshold in the epoch 1 data set, an image of the residuals at the 
source location is included as part of Figure [SI for completeness sake. No observable position 
offsets were detected between the sources detected in epoch 1 and those detected in epoch 
2. 

Comparisons between the VLBl images and the higher frequency ATCA images were 
complicated by a number of factors. Firstly, source identiflcation in the ATCA images was 
limited by the large-scale diffuse emission and confusion as a result of the signiflcantly lower 
resolution of these images compared to that of the VLBl images. Secondly, when compared 
to the VLBl image, a systematic offset of 900 mas, 160 mas and 205 mas was observed in 
the 8.6 GHz, 17/19 GHz and 21/23 GHz images, respectively. The 4.8 GHz ATCA image 
was not used in any further comparisons in this paper as a result of its poor resolution. 
The large offset observed at 8.6 GHz results from a 945 mas error in the phase calibrator 
position of the original 1993 observation. The ~ 200 mas offset observed between the high 
frequency ATCA images and the VLBI images is well below the ATCA beam size and is 
typical given atmospheric conditions at these frequencies and the 5° separation between the 
phase calibrator and target source. 
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To aid in discrete source identification, four large elliptical Gaussian components were 
used to model and effectively subtract the large-scale diffuse emission from each of the ATCA 
images. The full-width half-maximum (FWHM) size, position angle (PA) and flux density 
of these components, as determined separately for each frequency observed with the ATCA, 
are listed in [31 To compare the resulting ATCA residual sources against the VLBI sources, 
the elliptical Gaussian components used in the VLBI model were transferred to each of the 
ATCA data sets and their positions adjusted accordingly to account for the systematic offsets 
previously mentioned. Each of the model components were significantly unresolved with 
respect to even the smallest ATCA beam. With the size, position angle and location fixed, 
the flux density of each model component was varied to fit the diffuse-emission-subtracted 
ATCA data. The modelled flux densities for each source and data set are listed in Table [31 
No significant residuals were observed in the ATCA image or visibilties after the large-scale 
structure and source components were subtracted from the image. 



4. Discussion 

4.1. Radio spectra and free- free absorption modelling of the compact sources 



Lenc &: Tingayl (l2006l ) and iTingayl (120041 ) showed that the spectra of compact radio 
sources in NGC 253 exhibited a sharp downturn at 1.4 GHz, and to a lesser degree at 2.3 
GHz, as a result of free-free absorption. To test a similar scenario quantitatively with the 
compact radio sources detected in NGC 4945, an analysis was performed using the data 
tabulated in Table [31 In this analysis, the second epoch VLBI data was used to model the 
low frequency end of the spectrum and the ATCA data was used to model the high frequency 
end of the spectrum. Archival 8.64 GHz data were also included to model the mid-r ange 
spectrum. Observations of compact rad io sources in M82 (lUlvestad fc Antonuccilll994j ) and 
NGC 253 (lUlvestad fc Antonuccil 119971 ) have shown that for existing sources, little or no 
change in flux is observed over periods of order a decade. Here we have assumed that if a 
source existed in both the archive 8.64 GHz data and our current ATCA or VLBI data then 
it has not significantly faded or brightened in the intervening years between the observations. 
Where no 8.64 GHz data exist we do not constrain our models to a lower limit based on the 
detection threshold at this frequency as we do not know if the source existed at that point 
in time. 



In our analysis we used the same models and techniques described in [Lenc fc Tingay 



( I2OO6I ). However, for sources that were detected in only two observing bands a simple power- 
law spectrum was used to model the source spectrum, whereas the more highly sampled 
sources were also modelled with a free-free absorbed-power-law spectrum and a self-absorbed 
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bremsstrahlung (free-free absorbed) spectrum. For the large-scale emission, lower limits were 
not applied at the low frequency end of the spectrum as this emission is totally resolved by 
our high-resolution VLBI observations. 

The spectrum for each source was tested against each of the models using a reduced 
criterion to determine the best fit. 5*0, tq and a from the model fits are listed in Table 
m In all cases, the sources were modelled best by either a simple power-law spectrum or 
a free-free absorbed-power-law spectrum. The model for each source is shown against the 
corresponding measured spectrum in Figure [71 Error bars of 10% are shown for the flux 
densities at 23 GHz, 21 GHz and 2.3 GHz, whereas slightly larger errors of 15%, 15% and 
20% are estimated for the 19 GHz, 17 GHz and 8.64 GHz measures, respectively. The larger 
errors associated with the lower frequency ATCA observations take into consideration the 
effects of confusion and the increased levels of diffuse emission. The data shown in Table 
m and Figure [7] show that at least 9 of the 13 compact radio sources modelled are free-free 
absorbed. 



Using the same criteria as lLenc fc Tingayl (120061 ). of the 13 compact sources for which 
spectral indices could be determined in Table IH three have flat intrinsic power law spec- 
tra (a > —0.4) indicative of H II regions dominated by thermal radio emission (T), the 10 
remaining sources have steep int rinsic spectra (S), as normally associated with supernova 
remnants (iMcDonald et al.ll2002l ). The proportion of synchrotron sources (77%) compared 
to therma l sources (23%) in NG C 4945 is similar to that found in other nearby starburst 
galaxies. (IMcDonald et al.ll2002l ) found 65% o f the compact radio s ources in M82 were syn- 
chrotron sources and 35% thermal, whereas (ILenc fc Tingayl 120061 ) found that 60% of the 
detections were synchrotron sources and 40% thermal (after correcting source 5.805—38.92 
with a spectral index of -1.2). The comparatively lower proportion of thermal sources de- 
tected in NGC 4945 compared to M82 and NGC 253 is believed to be a result of the lower 
resolution of the ATCA compared to the VLA. Thermal sources are more prominent at the 
high frequency end of the spectrum, however, even at 23 GHz the ATCA does not have 
sufficiently long baselines to easily resolve thermal sources from each other and from the 
large-scale diffuse emission. 

The red circles in Figures El^a) and [31(b) show how the free-free opacity varies across the 
field of the nuclear region of NGC 4945. Although there appears to be significantly greater 
free-free absorption toward the nuclear region of NGC 4945, the ionised medium appears to 
be clumpy. As we are viewing the galaxy almost edge-on this may be the result of seeing 
sources in front of the disk, embedded within the disk, or behind the disk, rather than being 
indicative of variation in the ionised medium. However, it is also possible that the absorption 
may be localised around the source itself, such as has been observed in young, dense, heavily- 
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20031 ) ■ 



embe dded clusters (iTurner et al.lll998l : iKobulnicky fc Johnsorull999l : iJohnson fc Kobulnicky 



Based on the model parameters, the flux density of each source at 1.4 GHz, 5*1.4, is 
calculated and listed in Table HI It is expected that up to eight sources may be detected 
at 1.4 GHz assuming a detection threshold of approximately 1 mJy beam~^. Follow-up 
observations of NGC 4945 have been proposed with the LBA to further constrain the free- 
free parameters of the sources at lower frequencies. Such an observation would also help to 
highlight if any of the source flux density measures have been adversely affected by resolution 
effects since the LBA beam at 1.4 GHz would be substantially larger than at 2.3 GHz. 



4.2. Source variation 



The mean of the individual source flux density ratios between the two VLBI epochs is 
1.06, while the median is 0.98. If only relatively strong sources are considered (S'2.3 > 6 mJy), 
to exclude those that are significantly affected by deconvolution errors and measurement 
uncertainties, we observe a mean flux density ratio of 1.00 and a median of 0.985 between 
epochs. This suggests that there are no significant systematic differences between the fiux- 
density scales or general source fading or brightening between the two epochs. None of 
the eight sources in this partially restricted sample exhibited flux variations above the 10% 
errors of our measures. If we take the median ratio as an upper limit on the median fading of 
supernova remnants in NGC 4945 then this would equate to a 1.5% reduction in flux density 
over the 1.9 yr baseline between the two VLBI epochs, or a fading of < 0.8% yr~^. Further 
observations, over a longer period of time, would be required to confirm whether this fading 
is indeed real. 



4.3. Comments on individual compact radio sources 

4.3.1. H II regions 

The sources 27.9681-00.044, 27.4722-04.705 and 27.3861-07.554 are fiat-spectrum 
sources that may be associated with H II regions. All three sources are comparatively weak 
at 2.3 GHz and are resolved in the low resolution VLBI images. Only one of the sources, 
27.3861-07.554, is detected in the high resolution VLBI images. The source 27.4722-04.705 
is detected only in the second, more sensitive, epoch image and is not believed to be a new 
source. 
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It should be noted that the classification of the three sources as H II regions is not 
entirely consistent with their brightness temperatures (~ 2 x 10^ K) at the low frequency 
end of the spectrum. Significant errors in the measured spectral index of these particular 
sources may have resulted in their misclassification. These errors arise from poor sampling 
in the spectral energy distribution (SED), scatter in the SED as a result source modelling 
errors and underestimates in the total flux density of extended sources in VLBI images. 
Another possibility is that the VLBI detections are of one or more supernova remnants 
embedded within H II regions. These possibilities could be tested with further observations 
at intermediate and lower frequencies by improving the modelling of the SED of these sources. 

None of the fiat-spectrum sources appear to be significantly free-free absorbed, even 
those that appear to be relatively close to the nuclear centre. Similar characteristics are 
also observed in NGC 253 where only two fiat-spectrum sources ar e detected at 2.3 GHz , 



both of which are only weakly free-free absorbed at this frequency (ILenc fc Tingayi l2006l ) . 
Any curvature due to free-free absorption is difficult to detect in these sources since the 
low frequency VLBI observations are only sensitive to high brightness temperature sources, 
furthermore, the absorption acts in the same direction as the intrinsic spectrum. 



4-3.2. Steep-spectrum sources 

A total of ten steep-spectrum sources have been identified in NGC 4945. Two of the 
synchrotron sources, 27.5291 — 04.631 and 27.5607 — 04.779, in combination with 27.5402 — 
04.720, resemble a jet-like source and may be associated with a relativistic jet near the nucleus 
of NGC 4945 (section § 14. 4p . Both of the synchrotron sources are significantly affected by 
free-free absorption and to a similar degree (tq = 17 — 18) suggesting that they are in close 
proximity of each other and deeply buried behind an ionised screen. 

Three of the synchrotron sources 27.6164-03.373, 27.5734-03.793 and 27.3659-06.900, 
have shell-like structures which may be associated with supernova remnants and have shell 
diameters of 110 mas, 100 mas and 60 mas, respectively. Assuming an average expansion 
velocity of v =10,000 km s^^ relative to the shell centre , similar to that measured for the 



supernova remnant 43.31-1-592 in M82 (jPedlar et al.lll999l ). this implies ages of ~ 100(10^/f 



yr, ~ 90(10^/f ) yr and ~ 50(10'^/f ) yr for these remnants, respectively. The source 27.4646 — 
05.174 also has a shell-like structure in the VLBI image and is assumed to be steep spectrum 
source since it is not detected at higher frequecencies. This source has a diameter of 70 mas 
and, if it is associated with a supernova remnant, has an age of ~ 60(10^/f ) yr. 



Two of the synchrotron sources 27.7631 — 02.880 and 27.2963 — 07.073 are comparatively 
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weak and are fully resolved in the high resolution VLBI images. The remaining sources, 
27.4646 - 06.555, 27.4949 - 05.063 and the brightest source, 27.2710-06.592, all are partially 
resolved and of similar size (30 x 25 mas) but exhibit no significant structure. These may 
be young remnants with ages of ~ 25(10^/f) yr. 

4.4. A relativistic jet at the nucleus of NGC 4945? 

An intriguing feature of our VLBI images of NGC 4945 is the jet-like morphology of 
27.5291 — 04.631, seen in Figures |4]( a) and ll](b) and, at higher resolution in Figures [6](g) , 
El^h) andini^i). The morphology of this source is unlike that of any of the other radio sources 
seen in our images of the NGC 4945 starburst. Rather than being a compact or somewhat 
resolved discreet source, 27.5291 — 04.631 has an appearance strongly suggestive of the core- 
jet morphologies seen in arcsecond and milliarscsecond scale radio images of powerful AGN - 
radio galaxies and quasars. 27.5291 — 04.631 consists of a bright and unresolved component, 
located at the north-east end of an elongated structure of approximately 250 mas angular 
length and < 80 mas angular width (5 pc and 1.5 pc, respectively, at the 3.82 Mpc distance 
of NGC 4945). The elongated structure consists of at least 4 regions of enhanced brightness 
when imaged at low resolution (Figures Hl^a) and|l](b)). In the high resolution images in 
Figures ini^g), MM) and[n](i), the elongated structure appears to be almost more resolved, 
with no particularly bright, compact components standing out, aside from the compact 
component at the north-eastern end of the structure. This indicates a near-uniform surface 
brightness distribution along the length of the elongated structure. The images at the two 
observing epochs are also highly consistent in the features detected. In particular, high 
resolution images from the data at both epochs show a highly linear elongated structure 
close to the compact north-east component, but a clear apparent bend in the structure at 
the south-west end of the structure. 

What is the nature of 27.5291 — 04.631? As noted above, the morphological similarity 
of the source to a typical AGN core-jet source is striking. We explore this possibility below. 
Another possibility is that 27.5291 — 04.631 is the result of a chance alignment along the line 
of sight of multiple supernova remnants, similar to the other remnants detected in our VLBI 
image. While a possibility, this interpretation does not appear overwhelmingly likely. For 
example, having detected ~ 10 remnants in an angular area of 25 sq. arcseconds (Figure [31 
and assuming a uniform distribution of remnants in this area, the chance that 3 remnants 
(probably the minimum required to explain the 27.5291—04.631 structure) occur in an area of 
~ 0.1 sq. arcseconds is ~ 0.01. Further, the high resolution observations in Figure [6] provide 
a further filter on the data, only detecting emission from the highest brightness temperature 
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sources in the field. Along with 27.6164 - 03.373 and 27.5734 - 03.793, 27.5291 - 04.631 
was one of only three high brightness temperature sources detected. If we revise the above 
calculation with a total of 5 high brightness temperature sources detected, and again assume 
a uniform distribution of these sources over the area of the starburst region, the requirement 
that 3 of these 5 lie together to explain the structure of 27.5291 — 04.631 has a probability 
of only ~ 0.001. Even assuming an increase in the density of remnants toward the centre of 
the starburst does not make the chance alignment argument a plausible one. 

The interpretation of 27.5291 — 04.631 as a jet implies that the compact component at 
the north-east end of the elongated structure would coincide with a compact, massive object 
that provides the engine for the jet, a massive black hole being the likely candidate. In AGN, 
the supermassive black holes coincide with the dynamical centroid of the host galaxy. This 
does not seem to be the case in NGC 4945, as can be seen from Figure [3l 

The one-sided appearance of the jet suggests that the jet may be relativistic. We measure 
the jet to c ounter iet surface brigh tness ratio and find R > 22. Assuming the continuous slab 
jet model (IBlandford et al.l 119771 ) . a constraint can be placed on the product pCos9, where 
j3 is the speed of the radio emitting jet material (as a fraction of the speed of light) and 9 is 
the angle that the motion of the jet material makes to our line-of-sight. In this calculation 
we assume optically thin emission and a jet spectral index of a = —0.7 where oc z/". We 
estimate f3Cos6 > 0.52, using a value of i? > 22. Since we only have a limit on R, it is not 
possible to derive hard constraints on (3, given 6, or vice versa. A value of (3 = 0.9 implies 
6 < 55° whereas (3 = 0.52 implies 6 = 0° (i.e. perfectly aligned with our line of sight). 

An alternative explanation of the one-sided appearance of the assumed jet is that the 
counterjet is heavily free-free absorbed. Our modelling of the jet SED suggests that it is 
significantly affected by free-free absorption (r = 17— 18), the counterjet would presumably 
be affected to an even greater degree. In our images, the extent of the jet feature (~ 5 



19951) and NGC 1068 flGallimore et al.lll996 



(Peck et al. 


2003), NGC 4151 ( 


Mundell et al. 


Gallimore. Baum & O'Dea 


1996) which have 



jet features 9 pc, 50 pc, 260 pc and 1.3 kpc in extent, respectively. It is thus feasible 
that the counterjet is completely obscured by free-free absorption behind a dense ionised 
torus. Evidence for fr ee-free absorption in gaseous disks has been found in other Seyferts 
( iGallimore et al.lll999l ) and would be consistent with the observations in NGC 4945. 



If the jet interpretation is correct then this would be the nearest detection of such a 
feature to date. No such structures are see n in the two pro totypical starburst ga laxies, M82 
dPedlar et allll999l : iMcDonald et al.ll2002h and NGC 253 Jlenc fc Tin gavl l2006h. The next 



nearest detection of a Seyfert jet is in NGC 4151 at a distance of 13.3 Mpc (IMundell et al. 
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19951) 



This presents the opportunity to study jet physics at the highest possible spatial 



resolution. 

The jet interpretation for 27.5291-04.631 is easily testable with further observations. If 
monitored every 1—2 years, the motion of components away from the core may be detectable. 
In particular, two bright components appear separated from the core by approximately 30 
mas, at both epochs shown in Figure El We have measured the component separations from 
the core for these components and can put a limit on their motion, relative to the core, of 
~ 10 mas yr~^. For an intrinsic motion of (3 = 0.9 and an angle to the line of sight of 45°, we 
calculate an apparent motion of ~ 25 mas yr~^. Given this ballpark calculation, several years 
of monitoring would be required to confidently detect any motion of jet components. We 
are continuing to monitor NGC 4945 in order to test the jet hypothesis by directly detecting 
the transport of energy from the core along the jet. 



4.5. The supernova rate in NGC 4945 

A lower limit on the supernova rate for type II supernovae associated with massive stars 
can be estimated, based on the number of detected remnants, their size and an assumed 
expansion rate. We detect of order 10 compact radio sources with sizes limits of 0.5 — 2 
pc, based on a distance of 3.82 Mpc to NGC 4945, that exhibit steep intrinsic spectra 
often associated with synchrotron emission from type II supernova remnants. Assuming an 
expansion rate of 10,000 km s~^ (rate of expansion of FWHM of a Gaussian model), we 
estimate a supernova rate of > 0.1 — 0.4 yr~^. This is a lower limit because sources may 
have been missed as a result of being intrinsically weak, significantly free-free absorbed, or 
fully resolved. As our VLBI observations of NGC 4945 currently span a period of less than 
two years, no significant evolution is detected in the supernova remnants and their mean 
expansion rate remains unclear. It is therefore appropriate to rewrite the supernova rate 
limit as usn > 0.1(i;/10^) yr~^, where v is the shell radial expansion velocity in km s~^. 

It is interesting to note that three times as many VLBI sources are detected in NGC 
4945 compared to NGC 253. Yet it is assumed, by considering the effects of confusion at 8.3 
GHz in VLA observations, th at there may be of order 100 compact radio sources in NGC 253 



(jUlvestad fc Antonuccilll994l ). approximately 10 times the number detected in NGC 4945 at 
a similar frequency with the ATCA. In section § 14.11 we discussed the detrimental effects 
of confusion and diffuse emission on the thermal source counts in NGC 4945. The same 
problematic effects limit counts of compact steep spectrum sources with the ATCA but even 
more so as these sources will be significantly weaker at higher frequencies. So it is likely that 
there are a great deal more sources in NGC 4945 than our observations lead us to believe. 
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How this will effect the lower limit of the supernova rate is unclear and will depend on both 
the age and number of sources in this population. 

An upper limit on the type II supernova ra t e may be estimated by statistically mod- 
elling the effect of non-detections. iLenc fc Tingayl (120061 ) developed a Monte Carlo simulation 
based on a hypothetical population of supernova re mnants with a uni f ormly distributed peak 



lumin osity between 5 and 20 times that of Cas A (IWeiler et al.lll989l : lUlvestad fc Antonucci 



19911 ). This model has been adapted for NGC 4945 by scaling the flux density of the super- 



nova remnants to the distance of this galaxy, 3.82 Mpc. Similarly, simulation parameters 
were set accordingly for the two epochs of VLBI observations and are summarised in Table 
|5l Our simulations suggest that our VLBI observations should have detected all of the hy- 
pothetical supernovae that were simulated for between the two epochs assuming there were 
no significant absorption effects. Given this simple model we determine a supernova rate of 
< 1.62 yr^^ at a 95% confidence level. However, it is known that the nuclear region of NGC 
4945 is highly obscured and it has been shown, section § 14. H that the nuclear region is also 
situated behind a screen of ionized gas. In a second, more complete simulation, the effects 
of free-free absorption on the upper limit of the supernova rate are taken into consideration. 
For this test a median value of Tq = 11.8 was taken from the free- free modelling in § 14.11 
With this simulation only 11% of the hypothetical supernovae were detectable in the second 
epoch observation and the estimated supernova rate increased to < 15.3 yr~^. 

A further method of estimating the supernova rate is to estimate the median age of the 
supernova remnant population based on the observed flux variation between two epochs. In 
section § 14.21 we observed no significant variation in the detected sources above our error 
limits. However, if we assume that the median variation in the strongest of sources is real 
then we find an upper limit of ~ 0.8% yr^^ to the fading rate in NGC 4945 . If the supernovae 
are assumed to fade in time as the —0.7 power of time (jWeiler et al.lll986l ) then we estimate 
that the median age of the supernova remnant population of ~ 85 years. This would imply 
a supernova rate of ~ 0.12 yr~^, a value that falls within the limits determined above. 



4.6. The star formation rate in NGC 4945 



The star-formation rate (SFR) for stars with mass M > b Mc?^ of a star-forming galaxy 
can b e directly related to its radio luminosity L^, at wavelength i/ (jCondorull992l : iHaarsma et al 
20001 ) . This relation is understood purely from radio considerations by calculating the con- 
tribution of synchrotron radio emission from supernova remnants and of thermal emission 
from H II regions to the observed radio luminosity. For consi stency, we have taken t he same 
form of this relation and the same assumptions described in iLenc &: Tingayl (120061 ). Based 
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on the total integrated flux densities measured in the ATCA observations between 8.64 GHz 
and 23 GHz, we derive a SFR of 14.4±1.4 (Q/8.8) M© yr'^ for the inner ~300 pc of the 
nuclear region of NGC 4945, where the factor Q accounts for the mass of star s in the interval 
0.1 — IOOM0 and has a value of 8.8 if a Saltpeter IMF (7 = 2.5) is assumed (IHaarsma et al. 



2OOOI ). Similar measurements based on the 24 GHz continuum of NGC 253 give a comparable 
result of 10.4 ± 1.0 Mq yr'^ Jhenc fc Tineavlbooeh . 



Estimates of the SFR can be made based on the measured FIR luminosity /^f/j? since 
a large proportion of the bol ometric luminos ity of a galaxy is absorbed by interstellar dust 
and reemitted in thermal IR. ICondoru (119921 ) determined the relationship between the SFR 
of a galaxy and the FIR luminosity by modelling the total energy emission of a massive star 
and assuming that the contribution of dust heating by old stars was negligible. In NGC 
4945 this may be complicated by the presence of an AGN which, based on observations 
of the relat ively high HCN/CO ratio, may cont ribute as much as 50% of the bolometric 
luminosity (jSpoon et al.ll2000l : ICurran et al.ll200ll ). However, more recent observations of CO 
spectral lines in the nuclear region suggest th at star formation, rather than AGN activity, is 
the primary heating agent (IChou et al.l 120071 ) and so the SFR/FIR relationship should still 
provide a good estimate of the star formation rate. 

For NGC 4945, the far- infrared lumi nosity, LpjR, is 1.69 x 10^'^ Lc^ when adjusted for 
a distance of 3.82 Mpc (IRice et al.l Il988l ). Based on ICondonI (119921 ). this leads to a star 
formation rate of 1.5 Mq yr~^. Thi s is significantly low er than earlier estimates (8.3 Mq 



yr ^) based on total FIR luminosity (jPahlem et al.lll993l ) as a result of the revised distance 



to NGC 4945. However, the new estimate is comparable to the value 2 — 8 M© yr~^ derived 
from th e mass of ionise d gas inferred from radio recombination line observations from the 
galaxy (IRoy et al.ll2005l ). This estimate is also compatible with the lower limit on the star 
formation rate, > 0.19 Mq yr^^, derived from Ha emission (jStrickland et al.ll2004j ). 



For most g alaxies it is al so possible to directly relate th e star fo r matio n rate with the 
supernova rate (ICondorull992l ). Using the relation derived by ICondonI (119921 ) . the supernova 
rate limits of > 0.1{v/10^)iysN < 15.3 yr^^ (§113]) give the limits 2A{v/10^) < SFR{M > 
5Mq) < 370 Mq yr~^ on the star formation rate. While the upper limit is poorly constrained 
due to the limited number of epochs thus far observed, it is of the same order of magnitude as 
the nuclear SFR determined from FIR emissions (1.5 M© yr~^) and in agreement with that 
determined from radio luminosities alone (14.4 ± 1.4 {Q/8.8) M© yr~^). The higher radio 
estimate of the SFR, compared to the FIR-based estimate, may suggest that the Q factor 
has been over-estimated. A lower Q would suggest an increased proportion of massive star 
formation and would also be consistent with the high densities associated with the nuclear 
region. Nonetheless, it is encouraging that the independent estimates of SFR are only a 
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factor of a few different. 



4.7. Comparisons to the nearby starbursts M82 and NGC 253 



Witli tlie recent improvements in distance measurements, it appears that the distances 
to NGC 4945, M82 and NGC 253 are aU within a few percent of 3.9 Mpc. This is a dramatic 
change, by up to a factor of two, from previous estimates. Such a change in distance has 
an significant effect on luminosity estimates, which are proportional to the square of the 
distance, and to a lesser degree on spatial measures, which are linearly proportional to the 
distance. With the three galaxies now converging on a similar distance it is interesting that 
the star formation rates and supernova rates are also converging. Based on FIR luminosities, 
we find that the star formation rates of NGC 4945, NGC 253 an d M82 are now 1.5 Mg, 



YT'\ 1.8-2.8 Mq yr-i flLenc fc Tingavll2006h and ~ 2.0 Mq yr'^ flPedlar. Muxlow. fc Wills 



20031). respectively. Si milarly, the superno va rates are > O.lff /lO"^) yr ^, > 0.14(^/10'*) yr" 



fiLenc fc Tingavl l2006h and ~ 0.07 yr-^ flPedlar. Muxlow. fc Willsl 120031 ) . for NGC 4945, 
NGC 253 and M82, respectively, based on source counts and sizes. 



5. Summary 

We have imaged NGC 4945 at 2.3 GHz over two epochs using the LBA to produce 
the highest resolution images of the nuclear starburst region of this galaxy to date. The 
VLBI observations have been complemented with additional data between 8.64 and 23 GHz 
observed with the ATCA. We find the following results: 



• Fifteen compact radio sources were detected with the LBA and 13 were identified 
higher frequency ATCA images. 

• In the highest resolution image (15 mas beam), four supernova remnants are resolved 
into shell- like structures ranging in size between 60 and 110 mas (1.1 to 2.1 pc) in 
diameter. Assuming an average radial expansion velocity of v = 10,000 km s~^, the 
remnants are estimated to be approximately 25(10'^/i;) to 100(10^/f) years of age. 

• By combining flux density measurements from 2.3 GHz LBA and high frequency ATCA 
observations, the spectra for 13 compact radio sources were determined. The spectra 
of 9 sources were found to be consistent with a free-free absorbed power law and 4 with 
a simple power law spectrum. 
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• Ten of the 13 sources have steep intrinsic spectra normally associated with supernova 
remnants. The three remaining sources have flat intrinsic power law spectra {a > —0.4) 
indicative of H II regions but may have embedded supernova remnants owing to their 
high brightness temperature (~ 2 x 10^ K) at 2.3 GHz. 

• Based on the modelled free-free opacities of 9 sources, the morphology of the ionised 
medium in the central region of NGC 4945 is complex and clumpy in nature. 

• A type II supernova rate upper limit of 15.3 yr^^ in the inner 250 pc region of NGC 
4945 was derived from the absence of any new sources between epochs, taking into 
consideration the improved distance estimate for the galaxy, a median free- free opacity, 
and the sensitivity limits of two observations over a period of ~ 2 years. A type 
II supernova rate of > 0.1(t;/10^) yr~^ has been estimated based on an estimate of 
supernova remnant source counts, their sizes and their expansion rates. 

• A star formation rate of 2.4(w/10^) < SFR{M > 5Mq) < 370 M© yr'^ has been 
estimated directly from supernova rate limits for the inner 250 pc region of the galaxy. 

• The supernova rates and star formation rates determined for NGC 4945 are, to within 
a factor of two, similar to those observed in NGC 253 and M82. 

• A jet-like source is detected that is offset by --^ 1050 mas from the assumed location 
of the AGN based on H2O megamaser emission, the HNC cloud centroid, the K-band 
peak and the hard x-ray peak. 

Wc need to more stringently constrain the frcc-frcc parameters by observing a further 
epoch using VLBI at 1.4 GHz. Frequent observations will further constrain the upper limit 
on the supernova rate by increasing the proportion of supernova events that may be captured 
between epochs. Observations over a longer period of time (~ 10 years) will also further 
constrain the lower limit of the supernova rate by enabling the supernova shell expansion 
velocity to be measured, and allow us to search for motion in the suggested jet, thereby 
confirming its nature. 
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Table 1. Summary of observations. 



Observatory 


Frequency 


a 


5 


Date 


Config. 


Duration 


Bandwidth 


At 




MHz 


(J2000) 


(J2000) 






(h) 


(MHz) 


(sec) 


ATCA 


4800.0 


13'^5™28!00 


-49°28'12'.'00 


19/20 NOV 1993 


6A 


11 


128 


15 




8640.0 










11 


128 


15 


LBA* 


2252.0 


13''5™27!50 


-49°28'6y00 


13/14 MAY 2005 




9 


16 


2 




2268.0 










9 


16 


2 


ATCA 


17000.0 


13>'5™27?50 


-49°28'6'.'00 


18 MAR 2006 


EW367 


3 


128 


15 




19000.0 










3 


128 


15 




21000.0 










3.5 


128 


15 




23000.0 










3.5 


128 


15 


ATCA 


17000.0 


13'^5™27?50 


-49°28'6'.'00 


24/25 MAR 2006 


6C 


9 


128 


15 




19000.0 










9 


128 


15 




21000.0 










9 


128 


15 




23000.0 










9 


128 


15 


LBA^ 


2269.0 


13»'5™27!50 


-49°28'6'.'00 


21 MAR 2007 




11 


16 


2 




2285.0 










11 


16 


2 




2301.0 










11 


16 


2 




2317.0 










11 


16 


2 



^Australian LBA observation included 9 hours at ATCA and Mopra, Tidbinbilla, Hobart and Ceduna, and 6 hours at Parkes. 

''Australian LBA observation included 11 hours at ATCA and Mopra, Hobart and Ceduna, 10 hours at Parkes and 5 hours 
at Tidbinbilla. 
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Table 2. Map Statistics. 



Figure 


Source 


Frequency 


Synthesized Beam 


a 


Peak Flux 


Integrated Flux 






GHz 


(mas) 


(mjy beam^-*-) 


(mJy beam~^) 


(mJy) 


121a) 


J izo / — oU4o 


2.3 


13 X 10 


0.342 


520 


1100 




JlidY— oU4d 


2.3 


13 X 10 


0.292 


530 


1100 


m^J 




Z.O 


oU X oZ 


n 1 f;o 
U. 15Z 


OA 

Z4 


i5U 


|31,b) 


IN 4y40 


Z.O 


oU X oZ 


n HQ c; 
U.Uo5 


OQ 

Zo 


1 fin 


|4],a) 


z/.ozyl — U4.Doi 


2.3 


80 X 32 


0.152 




30 


|4],b) 


z r .ozyi — U4.Dol 


2.3 


80 X 32 


0.085 


6.8 


36 




Z(.OiD4 — Uv5. o/o 


Z.O 


oU X oZ 


n 1 c;o 
U. 15Z 


O 1 

z.x 


0.5 




Z(.u1d4 — (jo.oio 


Z.O 


Qfl \y QO 

oU X oZ 


n HQ c; 
U.Uo5 


1 Q 

l.cS 


Q O 

o.Z 




rsKjK^ 4y4o 


A Q 
4.0 


01 QO \y 1 (^Cn 

ZloZ X iDoU 


U.5ZU 


OoU 


ovnn 
Z ( UU 


Mo) 


InLtU 4y4o 


8.6 


1155 X y55 


0.202 


1 Tn 
i/U 


1300 




IN LjC' 4y40 


1 7 


fill \y c; 1 fl 
oil X 51U 


U.Z5U 


( 


y4u 


^^^^ 


rsKjK^ 4y4o 


1 n 

iy 


C;fiA A A'i 

5dU X 44o 


n ovn 
U.Z ( U 


OO 


c /I n 
o4U 


Me) 


InLtO 4y4o 


21 


495 X oyy 


0.320 


A O 

4o 


800 


Mi) 


iN LjC' 4y4o 


Zo 


/1c;q QfiQ 
45^5 X oDo 


n Q/1 n 
U.o4U 


AR 

40 


( oU 




Z(.uio4 — Uo.oro 


Z.O 


1 fi 1 c; 
ib X 15 


U. i4y 


1 Q 
i.O 


1 n 
lU 




zr.oio4 — Oo.olo 


2.3 


16 X 15 


0.110 


0.76 


9.0 




Z ( .Diu4 — {Jo.o 1 o 


Z.O 


1 fi 1 c; 
iD X 15 


U.U / 5 


1 n 
l.U 


n fi 

y.D 


M<i) 


0'7C:'7'3/1 HQ 'zriQ 

Z/.o ( o4 — Uv5. /yo 


2.3 


16 X 15 


0.149 


i.O 


11 


isle) 


OVCl'VJ/l HQ '70'? 


2.3 


16 X 15 


0.110 


2.2 


12 


Mi) 


Z r .0 i o4 — Uv). /yvS 


Z.O 


1 fi 1 c; 
iD X 15 


U.U / 5 


i.y 


1 o 
iZ 


IMs) 


Z/.ozyi — U4.Dol 


2.3 


16 X 15 


0.149 


O.O 


38 




Zf.oZyi — U4.0ol 


Z.O 


1 fi 1 t; 
iD X 15 


U. iiU 


o.U 


Qn 
v5U 




Z/.oZyi — U4.uol 


Z.O 


1 fi 1 c; 
iD X 15 


U.U / 5 


Q Q 
O.O 


Q /I 

v54 


Hj) 


27.4949 - 05.063 


2.3 


16 X 15 


0.149 


2.2 


5.7 


Hk) 


27.4949 - 05.063 


2.3 


16 X 15 


0.110 


2.4 


5.8 


Hi) 


27.4949 - 05.063 


2.3 


16 X 15 


0.075 


2.3 


5.7 


Elm) 


27.4646 - 06.555 


2.3 


16 X 15 


0.149 


2.5 


5.9 


En) 


27.4646 - 06.555 


2.3 


16 X 15 


0.110 


3.0 


7.8 


Eo) 


27.4646 - 06.555 


2.3 


16 X 15 


0.075 


2.8 


6.9 


Ep) 


27.4646 - 05.174 


2.3 


16 X 15 


0.149 


0.66 


2.3 


Eq) 


27.4646 - 05.174 


2.3 


16 X 15 


0.110 


0.82 


4.1 


Er) 


27.4646 - 05.174 


2.3 


16 X 15 


0.075 


0.70 


3.2 


Es) 


27.3861 - 07.554 


2.3 


16 X 15 


0.149 


1.1 


1.0 


Et) 


27.3861 - 07.554 


2.3 


16 X 15 


0.110 


0.91 


2.3 


Eu) 


27.3861 - 07.554 


2.3 


16 X 15 


0.075 


0.97 


1.6 


Ev) 


27.3659 - 06.900 


2.3 


16 X 15 


0.149 


1.1 


3.1 


Ew) 


27.3659 - 06.900 


2.3 


16 X 15 


0.110 


0.84 


5.7 


Ex) 


27.3659 - 06.900 


2.3 


16 X 15 


0.075 


0.93 


4.4 


Ey) 


27.2710 - 06.592 


2.3 


16 X 15 


0.149 


11 


37 


Ez) 


27.2710 - 06.592 


2.3 


16 X 15 


0.110 


9.1 


30 




27.2710 - 06.592 


2.3 


16 X 15 


0.075 


9.9 


33 
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Table 3. Source flux densities. 



Source 


FWHM Size^ 


P.A.b 


82.3(2005) 


32.3(2007) 


S8.64 


Sl7 


Sl9 


S21 


S23 




(mas) 


(degrees) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


27.9681 - 00.044 


155 X 43 


68 


6.8 


5.1 




3.4 


3.6 






27.7631 - 02.880 


91 X 44 


73 


1.6 


2.9 


17 


8.2 


6.8 


5.1 


4.0 


27.6164 - 03.373 


105 X 54 


44 


13 


12 


8.4 


6.0 


3.9 


3.0 


2.7 


27.5734 - 03.793 


87 X 56 


83 


18 


18 


13 


6.3 


7.2 


6.6 


4.9 


27.5607 - 04.779<= 


97 X 62 


45 


9.6 


8.6 


20 


8.1 


6.6 


5.6 


6.7 


27.5402 - 04.720'= 


156 X 12 


-55 


9.7 


15 












27.5291 - 04.631<= 


72 X 37 


73 


15 


14 




17 


16 


13 


11 


27.4949 - 05.063 


67 X 29 


72 


8.3 


8.2 


38 


20 


19 


14 


16 


27.4722 - 04.705 


85 X 37 


90 




2.9 


16 










27.4646-05.174 


77 X 62 


59 


9.1 


9.0 












27.4646 - 06.555 


79 X 33 


77 


11 


12 


15 


8.9 


9.0 


9.7 


8.8 


27.3861 - 07.554 


84 X 39 


-84 


3.4 


4.0 


2.8 


2.3 


1.3 


4.5 


6.1 


27.3659 - 06.900 


79 X 35 


72 


6.3 


6.4 


14 


7.5 


7.9 


7.6 


6.4 


27.2963 - 07.073 


78 X 39 


79 


5.2 


4.1 




1.3 


1.7 


2.4 


1.2 


27.2710 - 06.592 


54 X 24 


83 


36 


33 


34 


13 


12 


12 


10 


27.4833 - 04.570^1 


5764 X 982'' 


42 






470 


300 


280 


240 


230 


27.4602 - 05.845<i 


422 X 422'= 


-51 






17 


19 


16 


18 


18 


27.4525 - 05.245<^ 


566 X 311'= 


34 






130 


66 


59 


62 


57 


27.4371 - 05.695'' 


9485 X 3059'= 


41 






600 


430 


390 


390 


370 



^Gaussian component sizes measured against short-baseline 2.3 GHz LBA data except where noted otherwise. 

''Position angle of major axis measured east of north. 

'=27.5291 - 04.631, 27.5402 - 04.720 and 27.5607 - 04.779 are the apparent source, jet and termination of a jet like source. 

''Sources used to model diffuse emission observed in ATGA data. 

*= Sizes of large-scale diffuse emission sources measured against 17 GHz ATCA data. 
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Table 4. Parameters of tlie free-free absorption models for all detected sources. 



Source 


So 
(mJy) 


a 


To 


Type^ 


C b 
'-'1.4 

(mJy) 


27.9681 - 00.044 


5.99 


-0.19 




T 


5.62 


27.7631 - 02.880 


534 


-1.49 


22.9 


S 


0.004 


27.6164 — 03.373 


84.5 


-1.03 


6.16 


S 


2.87 


27.5734 - 03.793 


101 


-0.93 


5.57 


s 


4.74 


27.5607 - 04.779 


551 


-1.47 


16.9 


s 


0.081 


27.5402 - 04.720 












27.5291 - 04.631 


1110 


-1.45 


18.4 


s 


0.079 


27.4949 - 05.063 


608 


-1.19 


19.1 


s 


0.033 


27.4722 - 04.705 


0.99 


1.29 




T 


1.52 


27.4646 - 05.174 












27.4646 - 06.555 


67.1 


-0.66 


6.62 


s 


2.05 


27.3861 - 07.554 


3.31 


-0.02 




T 


3.33 


27.3659 - 06.900 


103 


-0.88 


11.8 


S 


0.23 


27.2963 - 07.073 


5.72 


-0.42 




S 


4.97 


27.2710 -06.592 


745 


-1.38 


11.2 


S 


1.84 


27.4833 - 04.570 


2150 


-0.70 






1690 


27.4602 - 05.845 


16.9 


0.02 






17.0 


27.4525 - 05.245 


875 


-0.89 






648 


27.4371 - 05.695 


1710 


-0.49 






1450 



^As discussed in the text. Thermal sources are shown as type T, and non-thermal syn- 
chrotron sources are shown as type S. 

''Estimated flux density at 1.4 GHz based on model parameters. 
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Table 5. The supernova rate upper limit based on Monte Carlo simulations run over two 

observing epochs. 



Epoch 


Time 


^obs 


Sensitivity 










(yr) 


(GHz) 


(mJy) 




(yr-^) 


(yr-^) 


1 




2.3 


0.91 








2 


1.9 


2.3 


0.51 


1.0 


< 1.62 


0.11 < 15.3 



*Monte Carlo test run with no modelling of free-free absorption. 

•^Monte Carlo test run with a median free-free opacity of Tq = 11.8 apphed. 

Note. — The supernova rate upper limit based on Monte Carlo simulations run over two 
observing epochs. The time between epochs, the observing frequency and sensitivity of the 
observation are listed. At the end of the second epoch the proportion of supernova remnants 
detected (Psn) in that epoch is listed together with an estimate of the supernova rate upper 
limit based on all observations prior to and including that epoch. 
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Fig. 1. — A three-col our composite, multi-wa velength view of NGC 4945. Red, green 
and blue indicate H a (IRossa fc Dettmarl 120031). K— band Two M icron All Sky Survey IR 



djarrett et aLll2003h. and Chand ra soft x-ray fISchurch et al.ll2002l ) respectively. Inset: HST 



Paa image ( IMarconi et al.ll2000h with ATCA 23 GHz contours overlaid, contours are loga- 
rithmic intervals of 2^/^, beginning at 1.02 mJy beam~^. 
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(a) J 1 237-5046 - 2005 epoch (b) J 1 237-5046 - 2007 epoch 




12^37^15,250' 15,245' 15,240' 15,235' 15,230' 12''37™15,250' 15,245' 15,240' 15,235' 15,230' 

Right Ascension (J3000) Right Ascension (J2000} 



Fig. 2. — Uniformly- weighted total-power map of the source used for phase calibration, 
J1237-5046, as observed with the Australian LBA at 2.3 GHz over two epochs. Map statistics 
for the individual maps are shown in Table [2l Contours are drawn at ±2°, ±22 , ±2\ ±22, • • • 
times the 3a rms noise. 
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(a) NGC 4945 - 2005 epoch 



. 37.9661-00.044 




a?, 7631-02. BBO 



27.6164-03 373 



27.5734-03.793 



27,5402-0- 
27.5607-04.7' 




-04.631, 



949^05.06'3 
7.4646-05.174 ■ 



,374646-06.55? 



17.3710-06. 593 
,36^-^)6,900 

s 27,2963-07,073 

27,3861-07,554 



27,6' 27,4' 
Right Ascension (J2000) 



(b) NGC 4945 - 2007 epoch 



28'00-- ^37,9661-00,044 




,37,6164-03 373 



27,5734-03,793 



27,5402-0. 
27,5607-04,7' 




-04.631, 
.27,47^2-04,705 
949^05,06'3 
7,4646-05,i74 



.<k 37, 3710-06, 593 

S9--06,900 

, 27.2963-07.073 



27.3861-07.554 



27.6' 27.4' 
Right Ascension (J2000) 



Fig. 3. — Uniformly-weighted total-power map of NGC 4945 as observed with the Australian 
LBA at 2.3 GHz over two epochs. Map statistics for the individual maps are shown in Table 
O Contours are drawn at ±2 2, ±2^, ±25, ■ ■ ■ times the 3a rms noise. The overlaid circles 
give an indication of the degree of absorption in the vicinity of the source, with the diameter 
of the circle being directly proportional to tq. The da shed blue circle and dashed black circle 
i ndicate the position al error of the Chandra x-ray (jSchurch et al.l |2002| ) and HST k-band 
(IMarconi et al.ll2000l ) peak, respectively. The solid blue circle a nd solid black circle mark the 
positio nal error of the B.9,0 megamase r (IGreenhill et al.l 119971 ) and the centre of the HNC 
cloud (ICunningham &: Whiteoakll2005l ). 



-34- 



(a) 27.5291-04.631 - 2005 epoch 



(b) 27.5291-04.631 - 2007 epoch 




O -04.B" - 




(c) 27.6164-03.373 - 2005 epoch 



(d) 27.6164-03.373 - 2007 epoch 



5 -03,40' 




13"05°37.63' 



Right Ascension (J2000) 



-49°38'03.25" - 



o -03,35" - 




13"05"27,63' 



37,63= 37,61 = 

Right Ascension (J2000) 



Fig. 4. — Uniformly-weighted total-power maps of 27.5291-04.631 and 27.6164-03.373 as 
observed with the Australian LBA at 2.3 GHz over two epochs. Map statistics for the 
individual maps are shown in Table O Contours are drawn at ±22, ±2^, ±2^, ■ ■ ■ times the 
3cr rms noise. 
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Ca) NGC 4945 - 4.8 GHz 



Cb) NGC 4945 - 8.6 GHz 



-49°27'55" - 




38.0^ 37.5^ 37. 

Right Ascension (J3000} 



(c) NGC 4945- 17 GHz 




13''05"'38.0'^ 37. 5"^ 27. 0=^ 

Right Ascension (J2000) 

(e) NGC 4945-21 GHz 



5 -28'05" - 




-49''37'55" - 




38. 0" 37.5^ 27. 0=^ 

Right Ascension (J2000> 



Cd) NGC 4945- 19 GHz 



-49°28'00" - 



O -gS'Oo" - 




13''05™38.0^ 37.5^ 37.0^ 

Right Ascension (J2000) 



CO NGC 4945 - 23 GHz 



-49°38'00" - 



13''05'"28.0^ 27. 5"^ 37.0" 

Right Ascension (J2000) 




IS'^OS^Be.O* 37.5^ 27 0^ 

Right Ascension (J2000) 



Fig. 5. — Uniformly-weighted total-power maps of NGC 4945 as observed with the ATCA at 
4.8 GHz, 8.6 GHz, 17 GHz, 19 GHz, 21 GHz and 23 GHz. Map statistics for the individual 
maps are shown in Table El Contours are drawn at ±2*^, ±2 2, ±2^^ ±22, • ■ • times the 3a rms 
noise for all maps except for (a) and (b) where the lowest contour is at 6cr rms noise. 
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Fig. 6. — Naturally- weighted total-power maps of compact sources detected in NGC 4945 
using all six antennas of the LBA at 2.3 GHz. Map statistics for the individual maps are 
shown in Table [2l Contours are drawn at ±2°, ±2^, ±2^, ±25, ■ ■ ■ times the 3a rms noise for 
all maps. 
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(p) 27.4646-05. 1 74 - 2005 epoch 



(q) 27.4646-05. 1 74 - 2007 epoch 



(r) 27.4646-05. 1 74 - Combined 




13'^5™57.475= 37.4n'0= 37.455^ 37.460' 

Right Ascension (J^OOO) 



1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 






% -05. 15" 






• 

■ , , , , 1 


, , , , 1 , , , , 1 , , , , ■ 


-05.30" 


■ , , , , 1 , , , , 1 , , , , 





37.470= 37.465' 37.450'' 

Right Ascension (J3000) 



170= 37.465' 37.450= 

Right Ascension (J3000) 



(s) 27.3861-07.554 - 2005 epoch 



Ct) 27.3861-07.554 - 2007 epoch 



(u) 27.3861-07.554 - Combined 



-49°3B'D7.50~ 






13''05"37 390' 37.366' 37.3B6' 37 364' 37.363' 
Right Ascension (J3000} 



13''O5"37.30O' 37.338' 37.366' 37.364' 37.3B3' 
Right Ascension (J3000) 



13''O5"37.30O' 37.3BS' 37.366' 37.364' S7.3B3= 
Right Ascension (J3000) 



(v) 27.3659-06.900 - 2005 epoch 



(w) 27.3659-06.900 - 2007 epoch 



(x) 27.365906.900 - Combined 



-49°3B'06.B5' 




Right Ascension fJ3000) 

Cy) 27.2710-06.592 - 2005 epoch 

I I I I ' ' ' I ' ' ' ' 

s> 





37.365' 37 360' 

Right Ascension (.13000} 



(z) 27.2710-06.592 - 2007 epoch 





(oa) 27.27 1 0-06.592 - Combined 
' ' ' ' I ' ' ' ' I ' ' ^ 

■0 CgiN 




13''05"3V.3V5= 37.370= 37.365= 
1 (J3000) 



13^05"37.375= 37.370' 37.365= 
Right Ascension [J3000) 



13^05"37.375= 37.370= 37.365= 
Right Ascension [J3000) 



6. — Continued. 
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Fig. 7. — Measured flux densities (symbols with error bars) and free-free absorption models 
(solid lines) for 17 detected sources (upper plots) and four diffuse sources detected with the 
ATCA data (lower plots). The error bars are 10% of the measured flux densities except at 
8.64 GHz where they are 20% of the measured flux density. 



